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In this Paper we report on our analysis of three Chandra observations of the accretion- 
powered millisecond X-ray pulsar IGR J00291+5934 obtained during the late stages 
of the 2004 outburst. We also report the serendipitous detection of the source in qui- 
escence by ROSAT during MJD 48830-48839 (July 26- August 4, 1992). The detected 
0.3-10 keV source count rates varied significantly between the Chandra observations 
from (7.2 ± 1.2) x 10~ 3 , (6.8 ± 0.9) x 10~ 3 , and (1.4 ± 0.1) x 10~ 2 counts per second 
for the first, second, and third Chandra observation, on MJD 53371.88 (Jan. 1, 2005), 
q '. 53383.99 (Jan. 13, 2005), and 53407.57 (Feb. 6, 2005), respectively. The count rate for 

?— i ' the third observation is 2.0±0.4 times as high as that of the average of the first two ob- 

\ servations. The unabsorbed 0.5-10 keV source flux for the best-fit power-law model to 

the source spectrum was (7.9±2.5) x 10~ 14 erg cm -2 s _1 , (7.3±2.0) x 10~ 14 erg cm~ 2 
s _1 , and (1.17 ±0.22) x 10~ 13 erg cm~ 2 s -1 for the first, second, and third Chandra 
observation, respectively. We find that this source flux is consistent with that found by 
ROSAT [w (5.4±2.4) x 10~ 14 erg cm~ 2 s" 1 ]. Under the assumption that the interstel- 
lar extinction, Nji, does not vary between the observations, we find that the blackbody 
temperature during the second Chandra observation is significantly higher than that 
during the first and third observation. Furthermore, the effective temperature of the 
neutron star derived from fitting an absorbed blackbody or neutron star atmosphere 
model to the data is rather high in comparison with many other neutron star soft 
X-ray transients in quiescence, even during the first and third observation. If we as- 
sume that the source quiescent luminosity is similar to that measured for two other 
accretion powered millisecond pulsars in quiescence, the distance to IGR J00291+5934 
is 2.6-3.6 kpc. 

Key words: stars: individual (IGR J00291+5934) — accretion: accretion discs — 
stars: binaries — stars: neutron — X-rays: binaries 
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1 INTRODUCTION 

The evolutionary link between the millisecond radio pul- 
sars and the low-mass X-ray binaries was established by 
the detection of the first accre tion-powered millisecond 
X-ra y pulsar SAX J1808.4-3658 iWiinands fc van der Klisl 
Il998l: IChakrabartv fc Morgan! Il99gft . Over the last years 
the number of known accretion-powered millisecond X- 
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ray pulsars has steadily increased. Recently, the discov- 
ery of the sixth member of this class was announced. 
The so urce IGR J00291+5934 was discovere d by INTE- 
GRAL jEckert et all 12004 IShaw et alJ 120051) on Dec. 2, 
2004 (MJD 53341). Using Ross i X-ray Timing Explorer 
(RXTE) follow-up observations, iMarkwardt et alJ <2004bft 
discovered pulsations at a frequency of 598.88 Hz. The pul- 
sar i s in a 147.4 min. orbit with a low-mass com panion 
star iMarkwardt et al-lEooial : iGallowav et al.ll2005lh A pre- 
viously undetected source was found at a magnitude of 
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R w 17.4 within the INTEGRAL error circle in an i?-band 
image obtained with th e Robotic Palomar 60-inch telescope 
l|Fox fc Kulkarnill2004h . Spectra obtained with the 4.2 m 
William Herschel Telescope of this variable source showed 
Ha and He II emission lines, securing the identification of 
this variable star as the count erpart of IGR J0029 1+ 5934 
jRoelofs et al.l 120041 : see also iFilippenko et alj 12004). Fi- 
nally, a variable radio and near infra-red sou rce was found 
at the position o f the optical counte rpart ijPoolevI Eool 
iFender et ai1l2004l : ISteeehs et al"ll2004l) . 

During outburst the X-ray, optical, and radio prop- 
erties of the accretion-powered X-ray pulsars are very 
similar to those of non-p ulsating neutron star soft X- 
ray transients (SXTs; e.g. IWiinands fc van der Klisl 1+9991 
for an X-ray account). In quiescence, the accretion pow- 
ered X-ray pulsars observed so far with Chandra and/or 
XMM-Newton have lower luminosities than those of many 
non-pulsating neutron star SXTs in quiescence. Further- 
more, in those quiescent systems where it has been 
possible to detect sufficient photons in X-rays to al- 
low for a spectral decomposition, the X-ray spectrum 
is significantly harder than for the non-pulsating neu- 
tron star SXTs 4U 1608-52, Aql X-l, SAX J1748.8-2021 
in the Globular Cluster NGC 6440, and XTE J 1709- 
267 in quiescence dAsai et al.lll996l: [Rutledee et al.ll2002al: 
Carrman^^^ljj2002i Jjonkerej L _aljj2003j: IWiinand s et alJ 
20051: lin't Zand et all l200ll : ICackett, et aljboOfl see also 
Jonker et alJ 12004a and references therein). However, 
there are systems that have i ntermediate s pectral prop- 
erties and luminosities (e.g. Ijonker et al.l 12004b') . Out 
of these sources with intermediate properti es Cen X- 
4 jRutledee et alJ 1200 ri . GRS 1741.9-2853 <Muno et alJ 
I2003D . and SAX J1810.8-2609 ijonker et aD l2004bl) have 
not been observed in outburst with the RXTE Proportional 
Counter Array. For that reason it is unclear whether pulsa- 
tions were present in outburst or not. 

The spectrum of several neutron star SXTs in qui- 
escence is found to be well-fit by a neutron star at- 
mosphere (NSA) model sometimes supplemented with a 
power-law component. Especially in sources with a quies- 
cent luminosity near 10 33 erg s -1 th e spectrum is dom i- 
nated by a strong t herm al component ijj o nker et al . 2004a). 
iBrown et all dl998ft and lColpi et alJ <200lft showed that the 
neutron star core is heated by pycnonuclear reactions dur- 
ing/immediately after accretion episodes. In quiescence the 
heated neutron star core cools in X-rays. For a given core 
temperature, from the assumed thermal properties of the 
neutron star envelope, the cooling core spectral energy dis- 
tribution is determined by the neutron star atmosphere. 
In theory, an NSA model fit to the quiescent thermal X- 
ray component provides the means to measure the mass 
and radius of the neutron star and hence constrain the 
equation of state of matter at supranuclear densities. In 
practice, n umbers typical for a canonical neutron star were 
found (e.g. iHeinke et all 120031. rendering support for this 
interpretation. However, there is an ongoing debate about 
whether the tem perature of the therma l component is vary- 
ing or not (cf. iRutled ee et al.| l2002al : ICamoana fc Stellal 
2003 ; ICampana et al l2004ll . Differences in temperatures ob- 
tained within one quiescent epoch may be explained by 
changes in the neutron star atmosp here due to ongoing low- 
level accretion iBrown et alj l20021. Similarly, differences in 



the observed neutron star temperature obtained before and 
after outburst activity can be explained by differences in the 
amount of resid ual unburned light elements in the neutron 
star atmosphere iBrown et alJ200*3) . Large changes on short 
timescales would render it unlikely that the soft/thermal 
component is due to cooling of the neutron star, limiting 
the applicability of the NSA model fit. Finally, there are 
currently two sources known which returned to quiescence 
after an accretion epoch lasting several years (i.e. KS 1731- 
260 and MXB 1659-298). In those cases it is thought that 
the observed thermal spectral comp onent is caused by the 
cooling of the neutro n star crust iWiinands et al 1 120021: 
iRutledee et~aHl2002ph . The observed changes in the X-ray 
spectral properties are also ascribed to the c ooling down of 
the neutron star crust iWiinands et alJl2004l . 

In this Paper we present the detailed analysis of three 
Chandra observations of the accretion powered X-ray pulsar 
IGR ,100291+5934 obtained during the late stages of the 
2004 outburst. 



2 OBSERVATIONS, ANALYSIS, AND 
RESULTS 

We observed IGR J00291+5934 three times after the 2004 
discovery outburst with the back-illuminated S3 CCD-chip 
of the Advanced CCD Imaging Spectrometer (ACIS) detec- 
tor on board the Chandra satellite. The observations started 
on MJD 53371.88 (Jan. 1, 2005; Chandra observation ID 
6179), MJD 53383.99 (Jan. 13, 2005; observation ID 6180), 
and MJD 53407.57 (Feb. 6, 2005; observation ID 6181; all 
times are in UTC). The net, on-source exposure times were 
4.7, 9.0, and 12.9 ksec, respectively. We limited the read- 
out area of the S3-chip to l/8th of its original size yielding 
a smaller exposure time per CCD frame in order to avoid 
pile-up. 

After the data were processed by the Chandra X-ray 
Center (ASCDS version 7.5.0), we analysed them using the 
CI AO 3.2.1 software developed by the Chandra X-ray Cen- 
ter. We searched the data for background flares but none 
were found, hence we used all data in our analysis. We ex- 
tracted data from a circular area with a radius of three pixels 
(~ 1.5") centred on the best-fit source position for the first 
two observations. This area encloses 90-95 per cent of the 
energy since, as we will show below, the source spectrum is 
not very hard. The reason why we use this small area is to 
exclude as much as possible background photons. In such a 
small area one expects 1±1 background photon for the first 
observation. During the third observation (MJD 53407.57) 
the count rate is higher than in the previous two observa- 
tion (see below), which led us to use a six pixel (~ 3") radius 
circular area centred on the source position for the source 
extraction. 

In the three observations we detect respectively 36, 63, 
and 190 photons from a position consiste nt with that found 
in the optica l bv | Fox fc Kulkarnl (l2004h and in X-rays by 
iNowak et aD J2004h . This yields 0.3-10 keV count rates of 
(7.2 + 1.2) x 10" 3 , (6.8 + 0.9) x 10 -3 , and (1.4 + 0.1) x 10~ 2 
counts per second. These background subtracted (see be- 
low) count rates are so low that pile-up is not a concern. 
The count rate for the third observation is 2.0 + 0.4 times as 
high as that of the average of the first two observations. Due 
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to the low number of source counts the Cash statistic was 
used for th e spectral analysis of the first two observations 
l)Cashll97flFl . This means that the background should not be 
subtracted when fitting the spectrum in order to maintain 
the counting statistics (although with a background as low 
as that in our case, this should not make much difference). 
For the extraction of the background spectrum in the third 
observation, we used an annulus with radius of 5-15". We 
rebinned the source spectrum of the third observation such 
that each bin contains at least 10 counts. We fitted the spec - 
tra in the 0.3-10 keV range using XSPEC <Arnaudlll99eh 
version 11.3.1. Since 10 counts per bin is still low for the use 
of x 2 statistics in the spectral fits of the third observation, 
we compared the results using the \ 2 statistics with those 
obtained using the Cash statistics; the results are the same 
within the errors. 

For each observation separately we searched for variabil- 
ity in the lightcurve. Only in the last observation (ID 6181) 
is there evidence that the source count rate was higher to- 
wards the end of the observation than near the start of the 
observation (see F igure[Q. A Kolmogorov-Smirnov test (see 
IPress et al.lll992T) testing the hypothesis that the count rate 
of the source is constant showed that the probability that 
the count rate is constant is 0.36, 0.36, and 8.4xl0 -5 for 
observation with ID 6179, 6180, 6181, respectively (i.e. the 
hypothesis is observed to be disproven at the 4.1<r level for 
the third observation). We searched for spectral variability 
within an observation by making hardness ratios for the first 
and second half of the total number of photons detected dur- 
ing each observation and by investigating the photon energy 
as a function of photon arrival time. The hardness is defined 
as the count rate ratio between the 0.3-1.5 keV band and 
the 1.5-10 keV band. There is no evidence for significant 
changes in the hardness ratio during any of the observa- 
tions. To check whether the hardness ratio changes between 
the observations we calculated the hardness ratio from all 
photons detected during one observation. The hardness ra- 
tio during the first, second and third observation is 2.8±1.1, 
1.3±0.3, and 2.3±0.4, respectively. The hardness ratio of the 
third observation is larger than that of the second observa- 
tion, albeit at the 2 a level. 

In order to check whether the spectrum changed be- 
tween the observations we fitted the spectra for each of 
the observations separately. We fit the spectra using sev- 
eral functions that are often used for quiescent neutron 
stars: a power-law, an NSA, and a black-body model, all 
of these were modified by the effects of interstellar absorp- 
tion. In these fits the value of the interstellar absorption, 
Nh, was fixed to 2.8 x 10 21 cm~ 2 similar to the value found by 
iNowak et alJ ll2004h . When we left the Nh as a free parame- 
ter in the spectral fits for the third observation the Nh was 
consistent with 2.8xl0 21 cm~ 2 . As an indication of the influ- 
ence that fixing the Nh has on the fit parameters we include 
the results of a blackbody, a power-law, and a blackbody fit 
to the first, second, and third observation, respectively, in 
Table 1. For the NSA model the neutron star radius and 
mass are fixed to 10 km and 1.4 Mq, respectively. Further- 
more, we assumed that the neutron star magnetic field is 
so low that it is unimportant for the NSA modelling, leav- 
ing only the normalisation (i.e. one divided by the source 
distance squared) and the effective temperature as free pa- 
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Figure 1. The 0.3—10 keV background subtracted lightcurve 
for the accretion— powered millisecond X-ray pulsar 
IGR J00291+5934 detected with the Chandra satellite on 
MJD 53407.57. The data were binned in 1000 second intervals. 
Variability in the count rate as a function of time within the 
observation is clearly detected. 

rameters in this model. The results of these fits are shown 
in Table 1. 

The likelihood ratio method, which is the basis of 
the Cash statistic, does not provide a direct test to the 
goodness-of-fit whic h would help one to distinguish between 
different models fe.g. lBevington fc Robinsonll992ll . In order 
to obtain a handle on the goodness-of-fit we used Monte 
Carlo simulations. They consist of simulating 10 4 counts 
spectra. Each spectrum is drawn from a parent distribution 
with Poisson statistics. The parent distribution is allowed 
to vary according to the covariances of the best-fit model 
parameters. The C fit-statistic for each of the Monte Carlo 
simulations is compared with that of the best-fit model. Af- 
ter the Monte Carlo simulations have been performed the 
fraction of simulations that gave a lower fit statistic than 
that of the best-fit model is given. This percentage is called 
the goodness. If the model provides a good description of 
the data the goodness should be close to 50 per cent. Both 
a very low goodness and a very high goodness indicate that 
the model does not represent the data accurately. As can be 
seen in Table 1, for many of the single component models 
for which Nh was held fixed at 2.8xl0 21 cm~ 2 the goodness 
is not close to the nominal 50 per cent. However, in none of 
the case we can exclude the model at high significance based 
on the goodness. 

The spectral fits suggest that the spectrum of the source 
is harder during the second observation than during the first 
and third observation, as can be seen from a comparison be- 
tween the best-fit power-law index or the blackbody/NSA 
temperature of the observations. Besides the single compo- 
nent models we tried fitting a model comprised of a black- 
body and a power law modified by the effects of interstellar 
absorption to the data in order to assess the contribution 
of a power law to the spectrum. However, due to the low 
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Figure 2. Top panel: The quiescent spectrum of 
IGR J00291+5934 detected with the Chandra satellite on 
MJD 53383.99 (Jan. 13, 2005). The best fitting power-law 
model is over plotted. The fit was done on the unbinned data 
using the C— statistic, but for clarity the data points have been 
rebinned in the plot. Bottom panel: The quiescent spectrum 
of IGR J00291+5934 detected with the Chandra satellite on 
MJD 53407.57 (Feb. 6, 2005). The best fitting blackbody model 
is over plotted. In both cases the Njy was fixed to 2.8xl0 21 

_2 

cm . 



number of detected photons during the first and second ob- 
servation the fit parameters were unconstrained. A fit using 
such a fit function resulted in meaningful constraints only 
for the third observation (and only when we kept the power- 
law index fixed to 2) . The results of this fit are also listed in 
Table 1. We also tried to fit an NSA plus power-law model 
to the spectrum of the third observation. However, with the 
power-law index fixed to a value of 2 the power-law did not 
contribute significantly to the fit (the normalisation was a 
factor of 10 3 less than that found for the blackbody plus 
power law fit). In Figure |5] we have plotted the spectrum 
showing a power-law fit to the dataset of the second Chan- 
dra observation and a blackbody fit to the third Chandra 
observation. 

The unabsorbed 0.5-10 keV source flux for the first, 
second, and third Chandra observation using the best-fit 
power-law model is (7.9 ±2.5) x 10~ 14 erg cm" 2 s" 1 , (7.3 ± 
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Figure 3. The 0.5-10 keV unabsorbed 2004/2005 outburst 
lightcurve for the accretion-powered millisecond X-ray pulsar 
IGR J00291+5934 as detected with the RXTE and Chandra satel- 
lites. The RXTE 2.5-25 keV observed flux (Galloway et al. 2005) 
is converted to an unabsorbed 0.5-10 keV flux using a power- 
law spectrum with index 2 and an interstellar absorption of 
N H = 2.8 x 10 21 cm" 2 . 

2.0) x 10" 14 erg cm" 2 s" 1 , and (1.17 ± 0.22) x 10" 13 erg 
cm" 2 s" 1 , respectively (see Figure^! the plotted uncertainty 
on these fluxes is at the 68 per cent confidence level). We 
also listed in Table 1 the unabsorbed 0.5-10 keV flux for 
each model for each observation. In between brackets we 
give the upper limit to the fractional contribution to the 
unabsorbed 0.5-10 keV flux for a power-law model in case 
of a blackbody or NSA model fit. We also give the upper 
limit to the fractional contribution for a blackbody and NSA 
model in case of a power-law fit. 

Finally, we extracted an 18 ksec. archival Positional 
Sensitive Proportional Counter (PSPC) ROSAT observa- 
tion of the Cataclysmic Variable RX J0028.8+5917 also 
known as DQ Her obtained over the period MJD 48830- 
48839 (July 26- August 4, 1992). Besides DQ Her, a source 
is detected at a position consistent within the R OS AT PSPC 
positional accuracy with that of IGR J00291+5934. There 
are 44 counts in a circular region with radius of 40" in the 
0.1-2 keV energy range. However, from a background deter- 
mination of a circular region 100" away, we find that 43 per 
cent of those counts can be attributed to the background. 
This yields a background subtracted source count rate of 
(1.4 ± 0.4) x 10" 3 counts per second in the PSPC instru- 
ment. If we use a power law with index 2.5 and fix the Nh 
to 2.8xl0 21 cm" 2 we get an unabsorbed 0.5-10 keV source 
flux of « (5.4 ± 2.4) x 10" 14 erg cm" 2 s" 1 at the time of the 
ROSAT observation. 



3 DISCUSSION 

We have presented the analysis of three Chandra obser- 
vations of the accretion-powered millisecond X-ray pulsar 



Table 1: Best fit parameters of the quiescent spectrum of IGR J00291+5934. NSA stands for neutron star atmosphere, BB refers to blackbody, and PL to power law. All 
quoted errors are at the 90 per cent confidence level. The value in between brackets in the flux column denotes the (95 per cent upper limit to the) fractional contribution 
of a PL in case of a BB or NSA model fit or a BB/NSA in case of a PL model fit. In case of an upper limit to the fractional PL contribution the PL index was fixed at 2. In 
case of an upper limit to the fractional BB/NSA contribution the BB/NSA temperature was fixed at 0.2 keV/10 6 K. For the first two observations the goodness-of-fit is 
expressed via the goodness parameter, whereas for the third observation the goodness-of-fit is given by means of the reduced \ f° r a certain degrees of freedom (d.o.f.). 



Obs ID 


MJD 
(UTC) 


Model 


N H X 10 21 
cm - 2 


BB radius 
(TMT ) 2 km 


NSA norm. 
, 10 kpc ,2 


PL norm. X 10~ 5 
phot. kcV- 1 cm- 2 3 -l 


Temp. / PL index 
BB (koV), NSA (log K) 


Unabs. 0.5- 
(crg cm - 2 


10 kcV flux 


Goodncas / 
per cent / d.o.f. 






v 10 kpc 1 


( d ' 


6179 


53371.88 


BB 


2.8° 


1 1+ 19 

-0.7 






0.27±0.05 


5.2X10- 14 


(< 43% b ) 


6% 


6179 
6179 
6179 


53371.88 
53371.88 
53371.88 


NSA 
PL 
BB 


2.8 a 
2.8 a 

°- a - 3 o.s 


°- 36 io.26 


0.21+ 1 - 01 
-0.17 


2.6 ± 0.8 


6.25±0.15 

3.2±0.7 
0.32±0.10 


5.5X10 -14 
7.9xl0~ 14 
3.5xl0~ 14 


(< 39% b ) 
(< 23/26% b ) 
(< 47% b ) 


34% 
3% 

31% 


6180 


53383.99 


BB 


2.8 a 


(5.8 + 4-S) X 10- 2 






0.54±0.09 


5.3X10 -14 


(< 68% b ) 


99% 


6180 
6180 
6180 


53383.99 
53383.99 
53383.99 


NSA 
PL 
PL 


2.8 a 
2.8 a 

0.7+ 2 - 4 
-0.7 




(0.40+° «) X lO" 2 


1.8 ± 0.5 

1.0+ 1 - 

-0.4 


6.65 + 0.10 
2.3±0.4 
1.7±0.6 


5.7xl0~ 14 
7.3xl0~ 14 
6.3xl0~ 14 


(< 53% b ) 
(< 15/18% i ') 
(< 10/10% b ) 


88% 
14% 
24% 



6181 
6181 
6181 
6181 



53407.57 
53407.57 
53407.57 
53407.57 



6181 53407.57 

a Parameter fixed. 
b 



BB 
NSA 
PL 
BB + PL 
BB 



2.8" 
2.8 a 
2.8 a 
2.8 a 
, + 2.7 
-1.3 



4-0.14 
5 -0.1 



1.16 + 0.57 
0.48^ 



8+1-5 
-0.3 



4.1 ± 0.6 

0.7+ 15 
-0.7 



0.31 + 0.03 

2.9 + 0.3 
0.28+0.05 / 2 a 
0.34 + 0.07 



Percentage contributed to the total flux / upper limit to the fractional BB/NSA or PL contribution to the total flux. 



(< 27% b ) 1.0/13 

9.0xl0~ 14 (< 25% b ) 0.96/13 

1.3xl0~ 13 (< 12/12% b ) 1.4/13 

1.0xl0~ 13 (67+33)% b 1.0/12 

6.5xl0~ 14 (< 25% b ) 1.0/12 
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IGR J00291+5934. The observations were obtained right af- 
ter the end of the 2004 discovery outburst. The count rate 
for the third observation is 2.0 ± 0.4 times higher than that 
of the average of the first two observations. This shows that 
even in quiescence, there is evidence for variability (at the 
~5 a level). However, owing to spectral variability between 
the three observations, the unabsorbed 0.5-10 keV source 
flux of the Chandra observations did not change significantly. 
We also detected the source in quiescence in a ROSAT ob- 
servation obtained in 1992. The flux that we measured in the 
Chandra observations does not differ significantly from that 
in the 1992 ROSAT observation. From this we derive that 
the source had reached quiescence already during our first 
Chandra obse rvation obtained ~ 31 days after the discovery 
of the source iShaw et alj|2005l) . 

Under the assumption that the interstellar extinction, 
Nff, does not vary between the observations, we find from 
a blackbody/NSA fit to the spectrum of the second obser- 
vation that the blackbody/NSA temperature is significantly 
higher than that during the first and third observation. Fur- 
thermore, a single model fit using either an NSA or a black- 
body model with the Nh fixed t o the value derived during 
outburst (N H = 2.8 x 10 21 cm~ 2 : lNowak et alJliooil yields 
a rather high temperature when compared with other neu- 
tron star soft X-ray transients in quiescence even for the first 
and third observation (cf. TBiackBody=0.27±0.05, 0.31±0.03 
keV/ log (T N sa) = 6.25 ± 0.15, 6.33±g;Je for the first 
and third observation of IGR J00291+5934, respectively. 
Whereas lo g Tnsa = 5.84 - 6.14 f or Aql X-l [range due to 
variability] . iRutledee et"ai]l2002al : TBi a ckBo dv=0-186±0.005 
keV / log (T N sa) = 5.99±0.02 for Cen X-4. ICampana et alJ 
12004: T B iackBodv=0.14±0 .02 keV for SAX J1810.8- 
2609, iJonker et af] l2004bl Tm as j ii ead y=0-24±0.02 keV for 
XTE J1709-267. IJonker et alJl2004at log (T NS a) = 6.02 - 
6.24 for SAX J1748.8-2021 in t he Globular Cluste r 
NGC 6440 [range due to variabilitvl . ICackett et alJl2005ft . 
Even when we leave the N_h as a free parameter the best- 
fitting blackbody temperature is 0.32±0.06 keV for the third 
observation. Such a high temperature can be explained if the 
source spectrum contains a harder component besides a sin- 
gle (absorbed) thermal component. For the last and longest 
observation, the signal-to-noise ratio is high enough to test 
this by fitting the data with a model consisting of a black- 
body and a power law, modified by the effects of interstellar 
absorption. This provides a good fit with a slightly lower 
blackbody temperature of 0.28±0.05 keV when we fix the 
power-law index to 2. In this case the power law contributes 
33 per cent to the source luminosity. 

iRutledee et al.l i2002al) and ICampana et al] (I2004D 
showed that the quiescent spectra and flux of Aql X- 
1 and Cen X-4 also vary during and between observa- 
tions. However, in their interpretation these authors differ. 
ICampana fc StelTal i2003T) showed that with the current data 
for Aql X-l, it is not possible to distinguish between intrin- 
sic source variability in the soft blackbody component and 
correlated variability in the interstellar column density and 
the power-law photon index. Here we show, for the first 
time, that besides in these two canonical neutron star soft 
X-ray transients, spectral and count rate variability also 
occurs in the quiescent X-ray emission of an accretion pow- 
ered millisecond X-ray pulsar. Significant variability in the 
lightcurve was also found during the last observation. At the 



95 per cent confidence level there is evidence for variability 
during an observation of the accretion-powered mil lisecond 
X-ray pulsar XTE J0929-314 JWiinands et all2005ll . As ex- 
plained in the Introduction, the variability observed in the 
quiescent flux of the neutron star SXTs KS 1731-260 and 
MXB 1659-29 iWiinands et al.ll200llwiinands et al.ll2004h 
likely has a different origin than the variability observed 
here. Similarly, the variability of SAX J1748 .8-2021 in the 
Glob ular Cluster NGC 6440 in quiescence llCackett et al] 
120051) can be caused by the fact that in between the two 
quiescent observations the source experienced an outburst. 
The quiescent flux of GRS 1 741.9-2853 was fou nd to vary 
by more than a factor of 5 jMuno et al] 120031. However, 
this source has, so far, never been observed in outburst with 
RXTE and hence it is unclear whether this source contains 
an accretion-powered millisecond X-ray pulsar. 

The ratio between the peak and lowest quiescent un- 
absorbed 0.5-10 keV X-ray flux is ~ 2 x 10 4 . This was 
found after correcting the peak flux given in lGallowav et all 
(2005) for the different band-pass assuming the peak-flux 
spectr um is similar to the spectrum found bv lNowak et al] 
<2004h . This ratio is low for a soft X-ray transient. Unfor- 
tunately, the distance to IGR J00291+5934 is not well con- 
strained, so in principle we cannot say whether the quies- 
cent l uminosity is high or the outburst luminosity is low (or 
both) . iGallowav et al] i200Etl mention a minimum distance 
of 4 kpc but as those authors also state, the uncertainty in 
this estimate is large. If we assume that IGR J00291+5934 
has a 0.5-10 keV quiescent luminosity of ~ 5 — 10 x 10 31 
erg s _1 , similar to that of the accretion powered mil- 
liseconds pulsars SAX J1808.4-3658 and X TE J0929-314 
iCampana et al]l2002t IWiinands et al.ll2005l1 . the distance 
to IGR J00291+5934 would be «2. 6-3.6 kpc. In light of the 
uncertainty of this estimate the values derived by us and 
Galloway et al. (2005) are close together. However, an NSA 
model does not provide a satisfactory fit to the spectrum of 
our third Chandra observation if we fix the model normalisa- 
tion to that obtained for a distance of 3.5 kpc (x r ed ~ 3-7 for 
14 degrees of freedom). If we let the Nh free the Xred = 1-5 
for 13 degrees of freedom with N H = 0.65±g;^ x 10 21 cm~ 2 
for a NSA temperature of 94±4 eV. 

Since this source is relatively bright in X-rays in qui- 
escence, a future, longer XMIM-Newton or Chandra obser- 
vation should provide better constraints on the spectral pa- 
rameters and test for the time scale of variability. 
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